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Abstract 

We simplify the formalism of Polder and Van Hove [Phys.Rev.B 
4, 3303(1971)], which was developed to calculate the heat transfer 
between macroscopic and nanoscale bodies of arbitrary shape, disper- 
sive and adsorptive dielectric properties. In the non-retarded limit, 
at small distances between the bodies, the problem is reduced to the 
solution of an electrostatic problem. We apply the formalism to the 
study of the heat transfer between: (a) two parallel semi-infinite bod- 
ies, (b) a semi-infinite body and a spherical body, and (c) that two 
spherical bodies. We consider the dependence of the heat transfer on 
the temperature T, the shape and the separation d. We determine 
when retardation effects become important. 

1 Introduction 

It is well known that for bodies separated by d » dw ~ ch/ksT the 
radiative heat transfer between them is described by the Stefan- Bolzman 
law: 

Q0h 3 c 2 V 1 2J 
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where 7\ and T 2 are the temperatures of solid 1 and 2, respectively. In this 
limiting case the heat transfer is connected with traveling electromagnetic 
waves radiated by the bodies, and does not depend on the separation d. For 
d < dw the heat transfer increases by many order of magnitude, which can 
be explained by the existence of evanescent electromagnetic field that de- 
cay exponentially into the vacuum. At present time there is an increasing 
number of investigations of heat transfer due to evanescent waves in connec- 
tion with scanning tunneling microscopy and scanning thermal microscopy 
(STM) under ultrahigh vacuum conditions PQI21E1- STM can be used for 
local heating of the surface, resulting in local desorption or decomposition of 
molecular species, and this offer further possibilities for the STM to control 
local chemistry on a surface. 

A general formalism for evaluating the heat transfer between macroscopic 
bodies was proposed some years ago by Polder and Van Hove [1 . Their the- 
ory is based on the general theory of the fluctuating electromagnetic field 
developed by Rytov jl] and applied by Lifshitz [3] for studying the conser- 
vative part, and by Volokitin and Persson [U] for studying dissipative part of 
the van der Waals interaction. In the present paper we significantly simplify 
the formalism of Polder and Van Hove. We reduce the problem to the finding 
of the electromagnetic field created by a point dipole outside the bodies. The 
formalism of Polder and Van Hove requires the determination of electromag- 
netic field for all space and for all position of a point dipole, and requires the 
integration of the product of the component of the electric and magnetic field 
over the volumes of two bodies. Our formalism requires only the evaluation 
of a surface integral over one of the bodies and is simplified further in the 
non- retarded limit (small distances between bodies), where the calculation 
of the heat transfer is reduced to the problem of finding the electrostatic 
potential due to a point charge. We apply the formalism to the calculation 
of the heat transfer between: (a) two semi-infinite bodies, (b) a semi-infinite 
body and a spherical particle, and (c) two spherical particles. Problem (a) 
was considered by Polder and Van Hove, but our derivation is much simpler. 
The problem (b) was recently studied by Pendry j2]. However he only con- 
sidered large separation and neglected retardation effects. However for large 
separations retardation effects becomes very important. 
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2 Formalism 



Following Polder and Van Hove [T], to calculate the fluctuating electromag- 
netic field we use the general theory of Rytov (see Ref. j3j). This method 
is based on the introduction of a fluctuating current density in the Maxwell 
equations (just as, for example, the introduction a "random" force in the 
theory of Brownian motion of a particle). For a monochromatic field (time 
factor exp(— iut)) in a dielectric, nonmagnetic medium, these equations are: 

V x E = i-B (2) 

VxH = -i^D+^j' (3) 

c c 

where, following to Rytov, we have introduced a fluctuating current density 
associated with thermal and quantum fluctuations. E, D, H and B are the 
electric and the electric displacement field, the magnetic and the magnetic 
induction fields, respectively. For nonmagnetic medium B = H and D = eE, 
where e is the dielectric constant of the surrounded media. Accordingly to 
Rytov the average value of the product of components of j* for local optic 
case is given by formula 

(j/(r,^)jf (iW)) = (ji(*)jk*(* , )) u 8(u-u') 

G/(r)if(r')) =A(T,u)u; 2 lme(u;)8(T-T')5 ik , (4) 

n(u) = z ,, = (6) 

and for nonlocal local optic 

(j{(v)j f k * (r0) w = A(T, W )^ 2 Im^(r,r', W ), (7) 

where Eik{r, r', u) is a non-local dielectric constant. For simplicity, in the 
derivation we will assume local optics. However, the same derivation is valid 
also for the non-local optics case, and the final results is the same in the sense 
that in both cases the problem of the heat transfer between two bodies is 
reduced to the problem of finding the electromagnetic field outside the bodies. 
Compared to Polder and Van Hove, our treatment includes non-local effects, 
such as the anomalous skin effect. 
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In order to calculate the radiative energy transfer between the bodies we 
need the ensemble average of the Poyting vector 



(S(r)) w = ( C /8tt) (E(r)xB*(r)) u + c.c, 



(8) 



at suitable point r. From the Maxwell equations it follows that the electric 
field, produced by random current density j f , is given by 



E l (v) = Jdr'D ik (r,r',u)j f k {r') 



(9) 



where function D ik (r,r' ,u) obeys the equations 

ViV fc - 6 ik (V 2 - (io/c) 2 )e(r)] D kj (r, r', u) = (Airuji/ c 2 )5 %J 5(r - r') (10) 

V;V fc - 6 jk (V 2 - (u/c) 2 )e(r')] D ik (r, r', u;) = (Attuji/ 'c 2 )5 l3 5{r - r') (11) 
The functions D ik (r,r',u) have the following symmetry properties [7] 

D ik (r,r',u) = D ki (r',r,oj) (12) 

The Poyting vector can be expressed trough the average products of the 
components of electric field. Using equations (|1UI11|) we get 



(E i (r)E j W u = J dr"A(r")^ 2 Ime(r")A fc (r,r") J D; fc (r' ) r") 



2 

C - J dv"A(r" 



VIV: - 5 lk V' 2 ) D u (r, r") - ^6 ik 5(v - r") 



D* k (r'y) 



D ik (r,r") 



(v£V? - S lk V' 2 ) D*,(r, r") + -^6 jk 5(v - r"j 
= jj drM(r")Vf|^(r',r") (VID u (t,v") - V?A*(r,r")) 
-D lk {v,v") (VlD^v'y) - V?L>* fe (r',r")) } - AunA^ReD^r, r') 

( A ( T i) - A(T 2 )y 

2i 



J dsdD lk (r,r") (VfD* fc (r',r") - V^(r',r' 



-D* k (r', r") (Vi'Afc(r, r") - Vj^r, r")) - 47raM(T 2 )ReA;(r, r'),(13) 



where we transformed the volume integral over bodies 1 and 2 to a surface 
integral over body 1. Assume that the two points r and r' lie outside the 
bodies. Using that for r/r' 

A fc (r,r")V'^(r',r") = Vj^r, r")^(r', r")) , (14) 
and performing surface integral in (j!3|) gives 

<^(r)^(r% = {m) ~ 2 f iT2)y J dS? ■ {D ik (v,v")V'D* jk (v',v") 

- D* k (r', r")V'D ik (r, r")} - 4nujA(T 2 )ReD l3 (r, r') (15) 
Using the Maxwell equation (J2J) we can write the Poyting vector as 

(S> u = jf- {V (E(r) ■ E*(r')> - ((E(r) ■ V) E*(r')> - c.c.} r=r , (16) 

In the non-retarded limit the formalism can be simplified. In this case 
the electric field can be written as the gradient of an electrostatic potential, 
E(r) = — V0(r). Thus the total Poyting vector becomes 

(^l = ^/dS-{([ExBl) w + C . C .} 
= ^/dS-{<-[V0xB*]> w + a C .} 
= ^-|dS-{([0VxB*]> w + c.c.} 

= ^ / dS ■ V «0(r)0*(r')>. - cc.) r=r , (17) 
In the same approximation we can write 

A fc (r,r / ) = --V 4 V^D(r,r / ), 

where the function D(r,r') obeys the Poisson's equation 

AL)(r,r') = -47r5(r -r') (18) 
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Using the identities 



D ik (r,v") (Vr^(r',r")-V^(r',r" 
i-V.jv^ [D(v,v") (V!D* k (r>,r") - V^* z (r',r»)) 

-D(r,r») (V?V^* fc (r',r'') - Vf Z^(r',r" 



D(v,v") (V^ fc (r',r")-V^(r',r' 

-lv i V;.L>(r,r")V; , L>*(r',r" 



formula (fT^j) gives 



(0(r)0*(r')) t 



(^(r)S i (r'> <i , = V < V;^(r)^(rO> t 
A(T X ) - A(T 2 ) 



2i 



/dS£ IT(r',r")Vfi}(r,r'') 



(19) 



(20) 



-D(r, r")V"D*(r'> r") ^ - 47rA(T 2 )Im£>(r, r) 



(21) 



3 Heat transfer between two flat surfaces 

In this section we apply the general formalism, developed in the previous 
section, to the problem of the heat transfer between two flat surfaces. This 
problem was considered some years ago by Polder and Van Hove pQ , and the 
non- retarded limit was recently studied by Pendry j2]. The present deriva- 
tion is more compact and demonstrate the elegance of our general approach. 
We suppose that the half-space z < is filled by a medium (temperature 
Ti) with reflection factors i?i p (q, uj) and Ri s (q,u) for s- and p- polarized 
electromagnetic fields, respectively, and the half-space z > d is filled by a 
medium (temperature T 2 ) with reflection factors -R2p(q, and i^q, w), and 
the region between the solids, < z < d, is assumed to be vacuum. Let q 
be the component of wave vector k = (q, p) parallel to the surfaces and 
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Since the problem is translationaly invariant in the x = (x, y) plane, the 
Green's function -Djj(r, r') can be represented by the Fourier integral: 

Ai(r,r') = J -le^^Ai^,/^^), (23) 
Taking the x- axis along the vector q, equations ()10j) and become 

(p 2 + ^)^^') = -^^-^) (24) 
(V + z') - iq-^D zx (z, z') = ~^8{z - z>) (25) 

p 2 D zx (z, z') - iq^-D xx (z, z') = (26) 

p 2 D zz (z, z') - iq^-D xz (z, z> = -^5(z - z>) (27) 
oz cr 

p 2 D xz (z, z') + iq-^D xx (z, z') = (28) 

Since the equations for D xy and D zy are homogeneous, these components of 
the Green function must vanish. Thus solving the system of equations (J2U) - 
(l2~%j) is equivalent to solving altogether two equations: equation (}2~4"j) for D yy , 
and the equation for D xx which follows from equations (f2o"|) and (|2^|) : 

if + &) Dxx(Z ' Z>) = ~^T 5{Z ~ Z ' ] - (29) 
D xz , D zx and D zz for z ^ z' can be obtained as 

iq Q \q Q g 2 <9 2 

D xz (T T\ ~.D XX ^ D xz ~,D XX i D zz -— — -D xx (30) 

In the vacuum gap, < z < d, the solution of (|24jl has the form 

2-ku j, 
pc 2 



D w (z, /) = e ^ z - z '\ + %,e ipz + i^e"^ 2 (31) 



At the boundaries z = and z = d, the amplitude of the reflected wave is 
equal to the amplitude of the incident wave times to corresponding reflec- 
tion factor. The Green's function D yy is associated with the s-wave, which 



7 



satisfies the boundary conditions 

v y = R ls (w y -^e ipz '^ for z = (32) 
w y = R 2s e 2ipd (v y - j^-e~ ipz J for z = d (33) 
Using (EU-O 

n ( P, P„ p^ipd ( p \p{z-z') _j — p -\p{z-z') 



pc 



.J p*F\* * I I \ 

J R ls e ip ^ +y ) + R 2s e 2ipd e-' ip{z+z,) 



(34) 



A s J 

A s = 1 - e 2ipd R 2s Ris (35) 



Equation (J29J) for is similar to equation (|24j) for and taking into 
account that _D XX is associated with p- wave it can be obtained directly from 
()34|) by replacing R s ^ R p : 



p C \ 2 



5 ra = (jj -> R P ] (36) 

We note that in our approach the calculation of the reflection factors for s- 
and p- waves is considered as separated problem, which , if necessary, can be 
solved by taking into account non-local effects. For the local optic case the 
reflection factors are determined by the well known Fresnel formulas 

_ Ejp — 3j _ p — 3j . . 

EiP+ Si p + Si 

where e$ is the complex dielectric constant for body i, 



lu 2 



T £ < - ( 38 ) 



c 2 



For homogeneous (in xy-p\ane) bodies the second term in (fi7)|) for z- compo- 
nent for Poyting vector can be written in the form 

Vi (E k (r)E:(r')) r=T> = V fc (E k (r)E* z (v)) = V 2 (E,{t)e;(t)) 
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The last quantity in this expression is purely real, and does not give any 
contribution to z - component of the Poyting vector. Thus, from equations 
(fT5Jl and (fTTIj) we obtain 



c± [A( Tl ) - A(T 2 )} f d 2 q [ / 9 2 , 

{ yv dzdz' yy[ ' ] 



167ra; 

^ D ± D ( zz >) 
dz : dz' 



(27T) 



([j/ -> s]) + c.c. [ 



(39) 



Using (jMj) and in (|S9*|l we finally obtain 

o7T JO Jq<u)/c 



(1- 




2 )(1- 


R 2p (u) 


2 )(ni(w) -n 2 (o;)) 




1 


- e 2i P d R lp (<j)R 2p {u;) 


2 



X 



2\p\d 



X- 



/ dujuj I d 2 ge~ 

JO Jq>w/c 

ImRi Juj) ImRor, (u) 



l-e-W d R lp (u)R 2p (u) | 2 



(ni(o») - n 2 (o>)) 



where 



m(u) = ( e hui/kBTl - 1 



(40) 



(41) 



is the Bose-Einstein factor of solid 1 and similar for n 2 . 

Fig. la shows the heat transfer between two semi-infinite silver bodies 
separated by the distance d and at the temperatures T\ = 273K and T 2 = OK, 
respectively. The s- and p-wave contributions are shown separately, and the 
p-wave contribution has been calculated using non-local optics (the dashed 
line shows the result using local optics). It is remarkable how important 
the s-contribution is even for short distances. The detailed distance depen- 
dence of S z has been studied by Polder and Van Hove within the local optics 
approximation, and will not be repeated here. The nonlocal optics contribu- 
tion to (S z ) p , which is important only for d < I (where I is the electron mean 
free path in the bulk), is easy to calculate for free electron like metals. The 
non-local surface contribution to lmR p is given by [Hj 



(ImR 



■p) 

surf 



j, U Q 
Up Kf 
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where £ (q) depends on the electron density parameter r s but typically £ (0) 
1. Using this expression for ImR p in (J4U)) gives the (surface) contribution: 



S, 



ekf 



surf 



10, 



2 kld*ti 



Note from Fig. la that the local optics contribution to S pz depends nearly 
linearly on 1/d in the studied distance interval, and that this contribution is 
much smaller than the s-wave contribution. Both these observations differ 
from Ref. 0, where it is stated that the s contribution can be neglected 
for small distances and that the p-wave contribution (within local optics) is 
proportional to 1/d 2 for small distances. However, for very high- resistivity 
materials, the p-wave contribution becomes much more important, and a 
crossover to a 1 /(independence of S pz is observed at short separations d. 
This is illustrated in Fig. lb and lc, which have been calculated with the 
same parameters as in Fig. la, except that the electron mean free path 
has been reduced from I = 560 A (the electron mean free path for silver 
at room temperature) to 20 A (roughly the electron mean free path in lead 
at room temperature) (Fig. lb) and 3.4 A (of order the lattice constant, 
representing the minimal possible mean free path) (Fig. lc). Note that 
when I decreases, the p contribution to the heat transfer increases while the 
s contribution decreases. Since the mean free path cannot be much smaller 
than the lattice constant, the result in Fig. lc represent the largest possible 
p-wave contribution for normal metals. However, the p-wave contribution 
may be even larger for other materials, e.g., semimetals, with lower carrier 
concentration than in normal metals. This fact has already been pointed 
out by Pendry: the p-wave contribution for short distances is expected to be 
maximal when the function 



ImRp f» 



Im- 



e + 1 



Im 



2^ P + 

LUp \U! p LOpT j 



is maximal with respect to variations in 1/r. This gives: 

2k B T 



where we have used that typical frequencies u ~ kgT/h. Since the DC resis- 
tivity p = in /(u p r) we get (at room temperature) p ~ 2nh/k B T w 0.14 f2cm. 
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To illustrate this case, Fig. 2 shows the thermal flux as a function of the sep- 
aration d between the surfaces when (^/r)/A; B T = 120 and hup/ksT = 15.5. 
One body is at zero temperature and the other at T = 273K. 

Fig. 3 shows the thermal flux as a function of the resistivity of the 
solids. Again we assume that one body is at zero temperature and the 
other at T = 273K. The solid surfaces are separated by d — 10 A and 
tkVp » ksT. The heat flux for other separations can be obtained using 
scaling ~ 1/d 2 which holds for high- resistivity materials. Finally, we note 
that thin high-resistivity coatings can drastically increase the heat transfer 
between two solids. This is illustrated in Fig. 4 shows the heat flux when 
thin films (~ 10 A) of high resistivity material p = 0.14 ficm, are deposited 
on silver. One body is at zero temperature and the other at T = 273K. 
(a) and (b) shows the p and s-contributions, respectively. Also shown are 
the heat flux when the two bodies are made from silver, and from the high 
resistivity material. It is interesting to note that while the p-contribution to 
the heat flux for the coated surfaces is strongly influenced by the coating, 
the s-contribution is nearly unaffected. 



4 Local heating of a surface by an STM tip 

It was pointed by Pendry |2j the local heating of a surface by STM tip can 
be used for local modification of a surface if the heat transfer is sufficiently 
great. To investigate the power of a hot tip to heat a surface Pendry modeled 
the tip as a hot sphere of the same radius R as the tip. This is a common 
approximation when calculating tunneling current and the same arguments 
justify its use for calculating heat tunneling. Pendry considered the case 
R << d « d w ~ ch/k B T and the electrostatic limit. However at large 
distances retardation effects can be important. In fact, it will be shown 
below that the heat transfer between a sphere and surface in the asymptotic 
limit (large separation) can be obtained directly from formula (|40p. 

Consider distances d << dw ~ ch/kBT (at T = 300 K we have dw 
~ 10 5 A). In this case we can neglect by the first integral in ()40j) . put p ~ iq, 
and extend the integral to the whole q - plane. Using these approximations, 
the second integral in ()40|) can be written as 

fj roo roc 

S z = — duu (mM - n 2 {u)) / dqqe- 2qd 

7T Z JO JO 
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f lmR lp (u;)lmR 2p (u;) \ 

x i 1 1 - e-^H^M i 2 + b sJ r (42) 

Now, assume that the medium 2 is sufficiently rarefied and consist of parti- 
cles with the radius R << d, and the polarizability a{uS). Then e 2 — 1 - * 
47ra 2 n << 1, where n is the number of particles per unit volume. Thus, 
when n — > it is enough to include only the first- non-vanishing terms in the 
expansion of the integrand of Eq . (|42j) in powers of e 2 — 1- The heat transfer 
between one particle and a surface can be obtained as the ratio between the 
change of heat transfer after displacement of body 2 by small distance dz, 
and the number of particle in a slab with thickness dz: 



% poo roc 

S z = 2- dcucu (m(w) - n 2 (cu)) / dqq 2 e- 2qd 

IX JO JO 

x ^ 2lmR lp (uj)lma 2 (uj) + ^—^j Imi?i s (cj)ima 2 (co> j (43) 

To treat this expression we assume that | s(u) \» 1 holds for all relevant 
frequencies. In the limit d <\ e \~ 1 / 2 dw, where e is taken at the characteristic 
frequency ~ ksT/h, the reflection factor of p-wave becomes 

Im^-^. (44) 

The polarizability of the sphere is determined by 

£ 2 - 1 n3 T 3Ime 2 
a 2 = — -R'\ lma 2 « -. (45) 

^2 + * I £ 2 | 

We describe the sphere the same dielectric function as the substrate: 

= 1 - ? j., y (46) 
u{uj + l/r) 

where r is the Drude relaxation time and io p the plasma frequency. In this 
case the p-wave contribution becomes 
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This result is in an agreement with calculations of Pendry To evaluate 
the s-wave contribution in the limit d <\ e | -1 / 2 dw, we use the integral 

[°° dqe-^lm^^ » f°° dglm^ 
Jo iq + s Jo lq + s 

= ~Im 1 1 So I & /_2 /2 d* (e" z - e' 3 j | = | | s I cos(0), 

where s = s(q = 0) = (u/c)y/e, \ s l 1 / 2 ,^ = args . Thus, 

S sz « 2 ■ lO 2 ^/^- 9 / 2 /? 3 ^^ 11 / 2 ^^)- 1 / 2 . (48) 

From the comparison (|4"7j) and (|48|) the s-wave contribution exceeds the 
]9-wave contribution for d > (dwc/uj 2 T) l l 2 ; for typical metals at room 
temperature, TiUp/ksT ~ 10 3 and ~ 10 2 so that d > 10 2 A. For 
dw | £ |~ 1/,2 < <i < dw we obtain 

S p2 « lQkTh- 1 ' 2 d- 2 R z c- 1 Ty 2 {ulrf^. (49) 

For g?vk | £ | _1 ^ 2 < < we obtain the s-wave contribution 

S S2 « lQkTn- 7l2 d-' 2 R z c- 1 T^ 2 {uj 2 p T)^l 2 . (50) 

Assume now that the spherical particle is so close to the surface that we 
can neglect retardation effects. In this case the problem is reduced to finding 
of electrostatic potential created by point charge, located in vacuum. Using 
image theorem jU] the electrostatic potential can be written in the form 

(r ' } = | r - r' | " ei + 1 | r - F | + 

t ± crw (« - (-D-^|S^M) cos L J_ (51) 

where we have choosen the origin of the coordinate system at the center O 
of the spherical particle, and taken the polar axis along the line connecting 
O with the center of "image" sphere O', and assumed that the points r and 
r' have the polar coordinates (r, 6, 0) and (r', 6', <$') with respect 0\ F is the 
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"image" of r' ; rj = r — h = (r*, 9i, <f>) where h is the vector connecting the 
centers O and O'. h = 2(R + d) and P™(cos6) is the associated Legen- 
dre function (see Fig. 5). At any interior point of the sphere the resultant 
potential is 



oo oo 



D(r,r') =J2J2 A n (r')r n Pn ( cos e ) cosm(0 - 0') (r < R) (52) 



m=0 n=m 



We expand the potential (JoT|) in the spherical harmonics around the center 
O using the formulas [TUj 

1 _ ~ ~ r n {n-m)\ 
I „ _/ i ' ' ' j 



\- - i m =on= m r' n+1 (n + m)\ 
x P™ (cos B) P™ (cos 9')e m cos m(<j> - (j)') (r < r') (53) 

|r-f'| ion^i 1 J (J + ro)!(n + m)! 

xP™(cos6)Pl n (cos6')e m cosm((i)- (j)') (r,r' < h/2) (54) 

( _ ir+m ^M = £(-l)^x4r(^ j^(cos^) (r < fc) (55) 

where e = 1, e m = 2 for m ^ 0. Using ([53~ J) -([5l) |) we can rewrite (jBTj) in the 
form 



m=0 ra=m ; =m 

a+n+i^i-1 r n {l + n)\ 



p-( cos 9)P™(C0S 9') 

v 7 ei + 1 #+"+1 (Z + m)!(n + m)!y 
+OrO(^ + (-l)' + " +I |^^('+3)i',"*(oos9)} en ,co S m W -«(56) 

Across the surface of the sphere the Green function D(r, r') (as the func- 
tion of r), and its normal derivatives must be continuous. These boundary 
conditions lead to the equation 



l=m 



yn+l ^ n _|_ m y 



+ ( 1} + " + £l + 1 h^ 1 (I + m)\(n + m)\ ) P ™( cos ^ 
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+t^-r +1 (;;:)^'>} <-> 

A n = (1 ~ £2 ^ ^ +1 (58) 

Outside the sphere the function C™(r') satisfies to Laplace's equation and 
can be expanded as 

oo , m v 

W) = E ( + J ^ m (cos^') (59) 



l=rn 



Using (|59|) in (|5*7j) and taking into account that equation (jHTjl must be satis- 
fied at arbitrary r' we obtain equations for coefficients and b 7 ^ 



a m 



»/ - ^»E L ; J I *"*"*"' + t^tt "r L, m a « ( 6 °) 



l'=m 



(n + m)\ E\ + I \ hj \n + m / 

''£l + l,tlV /J l v (/' + m)!(n + m)! + + m, 



If we make the replacement 



<i = K { ^%Ki + Ku (62) 
(n + mjl 



then from (|6Uj) and we obtain 



= (63) 
Let us introduce the dimensionless coefficients x n i 

n ' \ J (n + m)!(Z + m)! n ' ni ~ 

(~^) n+m , 7 A w A;X n ;, (64) 

{n + my.[l + my. 

where 

7? - £i ~ l 
lP ~ Ei + V 
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£ = (R/h), X n = —A n /R 2n+1 . The coefficients x n i obey the equation 



*„, = („ + m y. + (e _t$[ m)l m («) 

The solution of (jHSj) can be found by iterations and has the form 

~ 2 ,, +1 , (n + l')\(l' + l)\ 



oo oo 

>f'_l_1 _o/"_ 



, 02 V V .a'+i^'+ix \ (n + W + Q! (V + l'W' + iy. (m) 
lP hhJ ^ W (/'-m)!(/' + m)!(r-m)!(/" + m)! + "- (66) 



Using dOHj) and (f6*T]) formula (|5fi|) can be significantly simplified 

1 r' 



oo OO / 1 I\ ^| 

D(r,r>) = £ £C7V) — T + - Pr(cos^) e m cosm(0-0') (67) 



m=0 Z=m 

Using (jfiTjl in (^Tj) we obtain 

i xt \ x*i /nv / n ^ ((n + m)!C n m (r)C(r') 

(0(r)0 (,)>„ = - {niH - n,M)EE( (B _ m) , fi2 „ +1 

ImA 1 

x- ^-cosm(0-0') I -47rA(T 1 )ImD(r,r / ) (68) 

I A n r J 

Using (|68|) in (JT7j) for the heat transfer between a sphere and a flat surface 
we obtain 

^ „oo oo oo oo 



T= ^ IAJ 

tt 2 ^ ^ nt4^ (n-m)!(/-m)! 

I D 12 \ 



n — m 



The above formalism gives, in principle, an exact solution of the problem in 
the non- retarded limit. However, for d « R extensive numerical calculation 
are necessary, because in this case the series converge slowly. The numerical 
results will be presented elsewhere. In the present publication we only present 
an approximate solution of the problem. 

Using image theorem for e » 1 and for the points r and r' close to 
the surface of the sphere, in the first approximation in the expansion of 
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the electrostatic potential in the sum of the potentials created by the image 
charges we can write the potential in the form 

D(r, r') = - Ri pr ^-r- - 

~ R 2p-, TT + R lp R 2p] (70) 

r — r'- r — r ■ 



where 



2R — r' 
r- = ; — (x,V,z) 



r' = (x', y', —h—z') and f • = (x-, y[, — h— z^). The value of the surface integral 
(|2*T]) does not change if we assume that the potential has the form ()7()j) in all 
space outside a sphere. Thus using the Green's theorem we can convert the 
surface integral to a volume integral over all space outside a sphere. This 
volume integral can be easily calculated using the fact that outside a sphere 
the potential D(r, r') obeys the Poisson's equation with the point charges 
located at r = r', r = F, and r = f \ . Performing the calculation gives: 

Im (0(r)0V)L = ^ KH - na(w)) Re ji^^-^ - 

R *lp R 2 P | ~, _ ~ | + RlpRlp | r / _ f . | ~~ R lp R *2p | ~7 _ ~ | | ( 71 ) 

Using (J71|) in (|17j) for the heat transfer between a sphere and a flat surface 
we obtain 

S = — &ujuj (ni(uj) — n 2 (uj))lmRi p (uj)lmR2p(uj) 
2n Jo 

Jo [(2£) - cos#] 30 a 

where a = UpT/Air. If we assume that in an accordance with (J47|) every the 
elementary volume of the sphere gives the contribution to the heat transfer 

Ik 4 

dS pz « — ^trr V 1 (7? - T 2 4 W, (73) 

P leO^ + rf-.Rcos^ 3 ^ 3 1 2 V 1 27 ' K ' 

then after integration of (fT3|) over the volume of the sphere we obtain the 
result of the same order magnitude as (J72J) . Because for <i < _R < g?vk | s | -1 / 2 
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in the accordance with (|48|) the s-contribution of small particle does not 
depend on the separation d formula (|48|) is valid also for small separation d. 
From the comparison (|T2"|) and (J4"H|) we get that for the sphere close to the 
surface s-contribution dominates for 

d > 10- 3 (d w c/uj p ) 3/2 R- 2 (uj p T)- 3/2 - 

For "normal" metal at room temperature d > 10 A. 

5 Heating of a particle by an STM tip 

Let us now consider the heat transfer between an STM tip, which we again 
model by a spherical particle with radius R 2 and the polarizability a 2 (oo), 
and a spherical particle with radius Ri and the polarizability ctru) located 
on a surface. We consider the case of large separation d » Ri,R 2 , and 
neglect by influence of the substrate on the heat transfer. At large distances 
the thermal electromagnetic field radiated by particle 1 can be considered 
as the radiative electromagnetic field of a fluctuating point dipole p^ with 
ensemble average 

(pfpf*) = AtTjlmaxMS* (74) 
The electric field of this point dipole is given by 

E = [3n(n ■ p/) - p'] (I - %) e ikr - A; 2 [n(n ■ p') - p'] — , (75) 

where k = u/c and where n is a unit vector along the axis connecting the 
two particle. The rate at which a particle 1 does work on a particle 2 is 
determined by 

P = A dujuj\ma 2 {uj) (EiE*)^ (76) 



(EiE*) w = A(Ti,u)lmai(u) ^ + J- + -jL j (77) 

After absorption in the particle 2 this work is converted into heat. In the 
same manner we can calculate the rate of cooling of particle 2 using the same 
formula by reciprocity. Thus the total heat transfer between the particles will 
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determined by 

h f°° 

S = 77. \o / duu[ni(u) - n 2 M]ImaiMlma:2(a>) 
( Ztt ) Jo 



/6 2k 2 2k A \ . , 



Using for the polarizabilities of particles expression (J45j) we obtain 



—[Ti — * T 2 ], (79) 

where g 7 ^ = hc/k B Ti. 

6 Relation between heat transfer and friction 

The heat transfer studied above is closely related to the frictional stress 
between bodies in relative motion, separated by a vacuum gap. In the last 
years this "vacuum" friction has attracted a great attention in the connection 
with the development of the scanning probe technique [H1IE1IHJII21IIH1II1]- In 
jH] we show that the frictional stress between the bodies having flat parallel 
surfaces separated by a distance d and moving with velocity V relative to 
each other is determined by formula which is very similar to (|4()|) : 



l07T J JO Jq<u/c 



(1- 




R Xp {u) 


2 )(1- 


R 2p (uj) 


2 ) 




1 - e 2i P d R lp {u)R 2p {u) 


2 



+^ r du [ d Ve- 2|p|d 

47T d JO Jq>w/c 

lm.Ri p (uj)lmR2p(uj) I dn{uj) 



dn{u) 
du 



| 1 - e- 2 \p\ d R lp (u)R 2p (u) | 2 V du J 
+ [p -> s] (80) 

For the distances d « dw this formula is reduced to: 

n 



h poo poo 

a = / du dqq 3 e- 2 i d 
2ix l Jo Jo 
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lmRi p (uj)lmR2p( y uj) ( dn(u) 



1 - e-WRiMRnJu) I 2 V du 



\p^s}\, (81) 



which is very similar to the corresponding heat-transfer formula (|42j) . For 
d < dw£~ l l 2 , where e(u) is taken at the characteristic frequency ~ k B T/h, 
the p- wave contribution is given by [S] 

, l hV 

G V » 0.3 (k B T/hujp) 2 jr, (82) 

and for <iiy£ _1 / 2 < d < dw- 

a ~ A J^^bT ^_ (g3) 

P 27T 2 cPdjy frdjp CUpX' 

For d < dwE -1 / 2 , the s-wave contribution a s becomes independent of d: 

a s « O.STr-^-^^^TV^y. (84) 

For <i > cfjy£ _1//2 we have: 

a s « 3.88tt- 2 d- 2 c- 2 k B TTuj 2 pV. (85) 

From the comparison (|8*2*Jl and (jSlj) we find that a s > cr p for d > c/u 2 t. 
For typical metal at room temperature this corresponds to d > 1 A. This 
is in drastic constant to the (conservative) van der Waals interaction, where 
the retardation effects become important only for d > c/u p jH]. Finally we 
carry out the transition to frictional stress between a particle with the radius 
R « d and semi-infinite body in (JSTj). To do this, we assume as in Sec. 4 
that the body 2 is sufficiently rarefied, i.e., that the difference £2 — 1 is small. 
Keeping only the first non- vanishing terms in the expansion of the integrand 
of (J81j) in powers of these difference, we get formula similar to (|43|) : 

i-2qd 




x <{ 2Imi?i p (u;)Ima 2 (cj) + ImiZis^Imo^^ } (86) 

In the limit d <\ e | _1//2 dw the reflection factor of p-wave is determined by 
()44j) and in (|86j) the p-wave contribution is reduced to the formula, which 
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was obtained by Tomassone and Widom ^3]. For the spherical particle in 
this limit we get 

3 hV (k B T\ 2 . - - 

and for | e | -1 / 2 d^y < d < d w we get 

, F « l*r-*?£ (^) 5/2 ar"V 2 -'c-'tf . (88) 
For d <| e | -1 / 2 <iv(/, cr s is independent of d: 

A 9/2 



a s » 3.3 • lO'TT-^tiVa^a^R 3 



(89) 



and for | £ | x / 2 d^ < d < dw we get: 



a s - 3.54vr- 3 / 2 ^ ^ a^ V^ 3 - (90) 

From the comparison (|8T|l and (}89| we get that cr s > a p for 0.1(cdw/c r i) 1 ^ 2 < 
d <\ e | -1 / 2 dw- For normal metal at room temperature 10 2 < d < 10 3 A. 
For | e | -1 / 2 dw < d < dw c s ~ o>- 

To estimate er for i? >> d we use the same approach as in Sec. 4. We 
define the frictional stress between the elementary volume dV and the semi- 
infinite body as da = (3<j/A7rR 3 )dV , where a is given by an expression for 
a spherical particle for R « d. After integration over the volume of the 
sphere for d « R <\ e \~ 1 / 2 dw we get: 

a p ^A.l^{^\- x WR- (91) 

Because in this limit a s does not depend on d, it is still determined by ([90)1 . 
From (|50|) and (jPTj) we get that cr s > a p for 

d > 2.4 ■ lO^cdiy^i" 1 ) 576 ^" 273 - 
For "normal" metal at room temperature and R ~ 10 3 A, d > 10 A. 
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Recently, Dorofeev et.al. ^2] have observed Brownian motion of a small 
metal particle connected by a spring to a holder, and located in ultrahigh 
vacuum in the vicinity of a gold surface. It was observed that the particle 
performed a stochastic oscillatory motion increased as the particle approach 
the gold surface. It was suggested that this increased damping is due to 
the coupling to the fluctuating electromagnetic field. From (J5UJ) and (J3TJ) 
we can estimate the damping constant j/m = a/mV. For d ~ 10 2 A,i? ~ 
10 3 A, m ~ 1CT 11 — 1CT 13 kg, and for "normal" metal at room temperature we 
get 7 P /m ~ 10~ 17 s^ 1 and 7 s /m ~ 10~ 13 s _1 . However in the experiment .12 
j/m ~ 10 2 s _1 . Thus the contribution of a fluctuating electromagnetic field 
to the damping constant can not explain the observed experimental date. 
This result is in an agreement with our earlier conclusion [TT] . 

The fluctuating electromagnetic field is an origin of the frictional drag, 
observed between parallel two-dimensional electron systems [HI]- I n the fric- 
tional drag experiments a current is drawn in the first layer, while the second 
layer is an open circuit. Thus no dc - current can flow in the second layer, 
but an induced electric field occur that opposes the "drag force" from the 
first layer. Recently we used the theory of a fluctuating electromagnetic field 
to calculate frictional drag force between two-dimensional electron systems 
[I?] . For frictional drag stress we found formula, which is very similar to 
(IHUj) . We found that for modulation-doped semiconductor quantum wells 
retardation effects are not important under typical experimental conditions, 
supporting earlier calculations where retardation effects always have been 
neglected [03 EI] ■ A striking new result we found is that for systems with 
high 2D-electron density, e.g., thin metallic films, retardation effects becomes 
crucial and in fact dominate the frictional shear stress a. 

7 Summary and conclusion 

We have developed a general formalism for the calculation of the radiative 
heat transfer between two bodies of arbitrary shape and properties. The 
formalism has been applied to three different cases, namely the heat transfer 
between (a) two flat surfaces, (b) a spherical particle and a flat surface, and 
(c) between two spherical particles. For two flat solid surfaces we have pre- 
sented numerical results for several cases of practical importance, namely for 
two "normal" (high conductivity) metals (silver), two (high resistivity) semi- 
metals and two silver metals coated by thin layers (10 A) of high resistivity 
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material. For high-resistivity metals the p-wave contribution dominates, but 
for "normal" (high conductivity) metals we found the remarkable result that 
the s-wave contribution dominates even for short separation between the 
solids. We have pointed out the close relation ship between the radiative 
heat transfer between two solids and the vacuum friction [..] which occur 
when one of the solids slide relative to the other solid. The formalism de- 
veloped in this paper may can be generalized to treat the vacuum friction 
between bodies with curved surfaces. 
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FIGURE CAPTIONS 

Fig. 1 (a) The heat transfer flux between two semi-infinite silver bodies, 
one at temperature Ti = 273 K and another at T 2 = K. (b) The same 
as (a) except that we have reduced the Drude electron relaxation time r for 
solid 1 from a value corresponding to a mean free path vpT = I = 560 A to 
20 A. (c) The same as (a) except that we have reduced / to 3.4 A. 

Fig. 2 The thermal flux as a function of the separation d between the 
surfaces. One body is at zero temperature and the other at T = 273f^. With 
(h/r)/k B T = 120 and hu p /k B T = 15.5. 

Fig. 3 The thermal flux as a function of the resistivity of the solids. The 
solid surfaces are separated by d = 10 A and huj p » k B T. The heat flux 
for other separations can be obtained using scaling ~ 1/d 2 which holds for 
high-resistivity materials. 

Fig. 4. The heat flux between two semi-infinite silver bodies coated with 
10 A high resistivity (p = 0.14 ficm) material. Also shown is the heat flux 
between two silver bodies, and two high-resistivity bodies. One body is at 
zero temperature and the other at T = 273/T. (a) and (b) shows the p and 
s-contributions, respectively. 

Fig. 5. Spherical particle (origin O) above a flat surface and its "image" 
(origin O'). 
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